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Abstract 

We study the chemical bonds of small palladium clusters Pd^ {n = 2-9) sat- 
urated by hydrogen atoms using electronic stress tensor. Our calculation includes 
bond orders which are recently proposed based on the stress tensor. It is shown that 
our bond orders can classify the different types of chemical bonds in those clusters. 
In particular, we discuss Pd-H bonds associated with the H atoms with high coor- 
dination numbers and the difference of H-H bonds in the different Pd clusters from 
viewpoint of the electronic stress tensor. The notion of "pseudo-spindle structure" 
is proposed as the region between two atoms where the largest eigenvalue of the 
electronic stress tensor is negative and corresponding eigenvectors forming a pattern 
which connects them. 



Wave function analysis; Theory of chemical bond; Stress tensor; hydrogenated Pd 



1 Introduction 



Studying structures and energetics of small palladium clusters is of great importance as 
the first step toward understanding their catalytic properties. In Ref. |1|, the structures 
and physical properties of small palladium clusters Pd^ (n = 2-15) and several larger 
clusters have been studied using density functional theory (DFT) calculation. They have 
investigated their isomeric structures extensively and found many energetically nearly de- 
generate isomers. In Ref. [2j, based on the lowest energy structures of Pd^ (n = 2-9) 
found in Ref. |1|, the role of small palladium clusters in catalyzing dissociative chemisorp- 
tion of molecular hydrogen has been studied by DFT calculation. The results include the 
structures of the Pd clusters under full hydrogen saturation. As for Pde cluster, they have 
reported detailed analysis of sequential H2 dissociative chemisorption starting from bare 
Pde cluster to Pd6Hi4 cluster. Their conclusion of this work is that the capacity of small 
Pd clusters to adsorb H atoms is substantially smaller on average than that of Pt clus- 
ters, indicating that Pd nanoparticles are less efficient than Pt nanoparticles in catalyzing 
dissociative chemisorption of H2 molecules. Although this may be not so industrially en- 
couraging result for the Pd clusters, the obtained structures have interesting features from 
the viewpoint of chemical bonds. Our paper is a follow-up study of these papers to learn 
more about nature of chemical bonds in Pd^ (n = 2-9) and their hydrogen-saturated 
versions using the electronic structures obtained by quantum chemical calculation. 

In our analysis of chemical bonds, we use the electronic stress tensor. This method is 
based on the Regional Density Functional Theory (RDFT) and Rigged Quantum Electro- 
dynamics (RQED) [^-10 and has been applied to several molecular systems 11-17 . Our 
method includes recently proposed bond orders fTT] which are defined using the electronic 
stress tensor. One of our purposes is to show their usefulness in the Pd clusters. As far as 
metallic clusters are concerned, our analysis had been only applied to Pt clusters 13 and 
AI4 cluster [17j so the present analysis can be useful basis for further research using our 
stress tensor based analysis. Special interest in these Pd clusters is that there seems to be 
H-H bonds within the clusters. Some of these H-H bonds are considered to form after hy- 
drogen molecules are dissociatively adsorbed to the Pd clusters [2] . It would be intriguing 
to investigate whether these H atoms are bonded from the viewpoint of electronic stress 
tensor and, if bonded, how the bonding nature differs from that of the free H2 molecule. 

This paper is organized as follows. In the next section, we brieffy explain our quantum 
chemical computation method. We also describe our analysis method based on the RDFT 
and the RQED, including the definition of our bond orders. In Sec. [3} we discuss our 
results. In Sec. |3.1[ we analyze the chemical bonds of the hydrogenated Pd clusters using 
our bond orders. In Sec. |3.2| we discuss the chemical bond using the stress tensor with 
special emphasis on the Pd-H bonds associated with H atom with high coordination 



number and the bonds between H atoms. In Sec. |3.3| we discuss a way to improve our 
bond order definition by integrating energy density over some area. We summarize our 
paper in Sec. [4} 
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2 Theory and calculation methods 



2.1 Ab initio electronic structure calculation 

We perform ab initio quantum chemical calculation for Pd clusters and their hydrides. In 
this work, calculations are performed by Gaussian03 program package [18j using density 
functional theory (DFT) with Perdew-Wang 1991 exchange and correlation functional 
(PW91) [l9). The 6-31G** basis set with polarization functions [20] - [22] has been used for 
hydrogen atoms and LanL2DZ effective core potential [23] for Pd atoms. Optimization 
was performed without imposing symmetry. 



2.2 RDFT analysis 



In the following section, we use quantities derived from the electronic stress tensor to 
analyze chemical bonds of bare and hydrogenated Pd clusters. This method based on 
RDFT and RQED [3-10 provides useful quantities to investigate chemical bonding such 
as new definition of bond order ^Tj-^lSj. We briefiy describe them below. (For other studies 
of quantum systems with the stress tensor in a slightly diflFerent context, see Refs 



25-36 



See also Refs. 37,38 for related discussion on energy density.) 

The basic quantity in this analysis is the electronic stress tensor density ^^{f) whose 
components are given by 



4m ^ 



dx^dx^ dx^ 



dx^ 



dx^ dx^ 



(1) 



where {fc, /} = {1, 2, 3}, m is the electron mass, '0i(r) is the ith natural orbital and Vi is 
its occupation number. 

By taking a trace of ^'^(r), we can define energy density of the quantum system at 
each point in space. The energy density e^{r) is given by 



(2) 



k=l 



We note that, by using the virial theorem, integration of (r) over whole space gives usual 
total energy E of the system: J e^{f)dr = E. 

Regional chemical potential fiR [3j is calculated approximately using £^{f) [llj . 



" dNn ^ n(r) ' 



(3) 
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where n(r) is the ordinary electron density at r. Since electrons tend to move from high 
Illr region to low /llr region, the distribution of /llr maps the chemical reactivity. 
Now, we define bond orders as 6:f(r) or /llji at "Lagrange point" 
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The Lagrange 

point Tl is the point where the tension density r'^(r) given by the divergence of the stress 
tensor 



4m 



^^ (O^-I ^-^^^^(0 



^,(r)-A^*(r) 



dx^ 
dx^ 



(4) 



vanishes. Namely, r {tl) = 0. r (r) is the expectation value of the tension density 

operator r'^(r), which cancels the Lorentz force density operator L(r) in the equation of 
motion for stationary state |7|. Therefore, we see that r^{f) expresses purely quantum 
mechanical effect and it has been proposed that this stationary point characterizes chemical 
bonding flT]. Then, our definitions of bond order are 



'tAB 



and 



£'rABirL)/nAB{rL) 
^fHH(^i)/"'HH(rL)' 



(5) 



(6) 



One should note normalization by the respective values of a H2 molecule calculated at the 
same level of theory (including method and basis set). 

We use Molecular Regional DFT (MRDFT) package |39| to compute these quantities 
introduced in this section. Some part of the visualization is done using PyMOL Molecular 



Viewer program 24 . 



3 Results and discussion 
3.1 Bond order analysis 

The optimized structures for hydrogenated Pd clusters Pd2H2, Pd3H2, Pd4H8, PdsHio, 
Pd6Hi4, PdyHie, PdsHie and Pd9H22 are shown in Fig. [T] In the figure, atoms are con- 



nected when the Lagrange point (Sec. 2.2) is found between them and the bond is colored 



to show the magnitude of the bond order (eq. ([5])). Fig. [2] shows the exactly the same 
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structures with the different bond order 6^ (eq. ([6])). These structures are obtained by 
re-optimizing the structures reported in Ref. ||2]. We performed optimization with multi- 
phcity of 1, 3 and 5 for each cluster and adopted the one with the lowest energy, which 
turned out to be singlet for all the clusters. We did not find much difference between the 
structures in Ref. p] and our re-optimized ones. For later use, we report that we performed 
similar procedure for bare Pd clusters Pd^ (n = 2-9) starting with the structures obtained 
in Ref. |1|. In this case triplets have the lowest energy for all the clusters. 

To show features of the chemical bonds in those clusters collectively and to exhibit 
usefulness of our bond order definitions, we plot bond orders against the bond length for 
all the bonds in the bare and hydrogenated Pd clusters in Fig. [3| In addition to our bond 
orders and 6^, we plot using conventional bond orders: the Wiberg bond index [40|, 



atom-atom overlap-weighted natural atomic orbital (NAO) bond order |41 , 42 and the 
Mayer's bond order |43j. It is apparent that our bond orders have better correlation than 
the other conventional bond orders. Therefore, we only show our bond orders for the 
following analysis. 

In Figs. [4]and[5| we re-plot the bond order v.s. bond length for and 6^ respectively, 
this time distinguishing between different types of bonding. The types we consider are 
Pd-Pd, terminal Pd-H, two-fold Pd-H, three-fold Pd-H, four-fold Pd-H and H-H. Here, 
two-fold Pd-H bond is associated with the H atom bridging two Pd atoms like Pd-H-Pd 
and, similarly, three-(four-)fold Pd-H bond with the H atom bonded to three (four) Pd 
atoms. Each type of bonding is plotted by different marks in Figs. [4] and [5] and number of 
each type for each cluster is shown in Table [TJ We note that some of the H atoms which 
we classified as having two-fold and three-fold Pd-H bonds have a bond with H atom 
in addition. In detail, H(12) and H(13) in PdsHio and H(10) and H(21) in PdyHie are 
bonded to each other in addition to two Pd atoms, which results in having three bonds 
from each H. Similarly, H(10) and H(ll) in Pd6Hi4 are bonded to each other in addition 
to three Pd atoms, making four bonds from each H. 

From Figs. [4] and [5| we see that bonds in the Pd clusters can be classified by the different 
slopes on the bond order v.s. bond length plane. There are a slope that corresponds to the 
Pd-Pd bonds and two slopes for Pd-H bonds. There are three outliers that correspond 
to H-H bonds. The fact that Pd-Pd bonds, whether they are in hydrogenated clusters or 
in bare clusters, are on a single slope indicates that the character of the bonding is not 
affected much by the hydrogenation. They on average become longer and weaker upon 
hydrogenation but the relation between the bond order and bond length is unchanged. 

Pd-H bonds may be classified in two groups. One has shorter bond length (< 1.9 A) 
and higher inclination to which terminal Pd-H bond and two-fold Pd-H belong. Some 
bonds in three-fold and four-fold Pd-H also belong to this group. Another group has 
longer bond length (> 1.9 A) and lower inclination and consists only of three-fold and 
four-fold Pd-H. For convenience, we call the former group "A" and the latter "B". Closer 
inspection of these H atoms with high-coordination numbers tells that there is no H atom 



^ We list all bond orders and enlarged structures with atom numbering in the supplementary materials 
(Fig. SI and Table SI). 
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which has only group B bonds. The bonds stem from H atoms consist of those of group 
A or mixture of A and B. See Table [2] for the detail. It is interesting that there are Pd-H 
bonds which are longer than some of the Pd-Pd bonds. Such long bonds are found in 



four-fold Pd-H bonds, which will be investigated more in Sec. 3.2, We also discuss in 



Sec. |3.2| that group A is characterized by a "spindle structure" and the group B by a 
"pseudo-spindle structure" . 

We now would like to give somewhat more quantitative analysis of these slopes in the 
bond order v.s. bond length relation by fitting the data points to linear and exponential 
curves. We fit data points which belong to the shorter Pd-H bonds (group A), the longer 
Pd-H bonds (group B) and the Pd-Pd bonds to functions in the forms y = ax + b and 
y = cexp(— rfx). Here, each group has 169, 13 and 134 data points, and in the fitting 
functions, y stands for or 6^ and x stands for the bond length. We summarize the results 
in Table |3| Roughly speaking, both functional forms give good fits, and the differences in 
the fitting parameters confirm the existence of two slopes for the Pd-H bonds (e.g. the 
inclination of the linear fits, the parameter a, is about 8 times larger for the group A 
than the group B). On closer look, the fits to b^ are better performed by the exponential 
function especially for Pd-H bonds. In contrast, the use of exponential form does not 
improve the fits to 6^ much, and in fact, the linear fits are slightly better for the Pd-H 
bond group A and the Pd-Pd bond group. 

As mentioned above, H-H bonds are shown in Figs. |4] and [5] as isolated points from the 
slopes of Pd-Pd and Pd-H. They cannot be put on a single slope and appear somewhat 
irregular. Since the hydrogenated Pd clusters here are formed by dissociative chemisorp- 
tion of H2 |[2j, these H-H are considered to be not trivial bonding of the H2 molecule. 
Namely, they are formed on or within the clusters and characteristic to the Pd clusters. 
They all have much longer and weaker bonds than the H-H bond in the H2 molecule. 
(Note that b^ and 6^ are unity for the H2 molecule by definition. The bond length of 
H2 in our computational method is 0.748252 A.) We will investigate these H-H bonds in 
more detail in Sec. 13.21 The existence of such H-H bonds are most notable difference from 



hydrogenated small Pt clusters which have been investigated in Refs. 13,44] 



3.2 Stress tensor analysis of chemical bond 

In the previous section, we have seen that there are several interesting bonding patterns 
found in the hydrogenated Pd clusters. We investigate these bonds in detail via the 
electronic stress tensor analysis. 

Before we discuss the Pd clusters, we show how the chemical bond of the hydrogen 
molecule is expressed by the electronic stress tensor (Eq. Q). In Fig. [6| on the left 
panel, we plot the largest eigenvalue of the stress tensor and corresponding eigenvector on 
the plane including the internuclear axis. On the right panel, we plot the tension vector 
(Eq. Q), which is normalized and whose norm is expressed by the color of the arrows. The 
sign of the largest eigenvalue tells whether electrons at a certain point in space feel tensile 
force (positive eigenvalue) or compressive force (negative eigenvalue) and the eigenvector 
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tells the direction of the force. We can see that the region with positive eigenvalue spreads 
between the H atoms, which corresponds to the formation of a covalent bond. In that 
region, the eigenvectors form a bundle of flow lines that connects the H nuclei. Such a 
region, called "spindle structure" [8j, is clearly seen in the panel. On the right panel, the 
vanishing point of tension vector "Lagrange point" (Sec. 2.2) is found at the midpoint of 
the internuclear axis, which is quite reasonable. 

We will now turn to the hydrogenated Pd clusters. We flrst look at the Pd-H-Pd 
bridging bond in Pd6Hi4 as shown in Fig. [7| We see the spindle structure between Pd 
and H just like the H2 molecule mentioned above, indicating the covalency of the Pd-H 
bond here. As for the Pd-Pd bond, although we found a Lagrange point and flow of the 
eigenvectors connecting Pd atoms, the eigenvalue in the region between Pd atoms has 
negative value. This indicates that the interaction between these two Pd atoms is weak 
and the bond between them is different from a covalent bond characterized by a spindle 
structure. We shall call this pattern a "pseudo-spindle structure". Namely, the pseudo- 
spindle structure has similar eigenvector flow to that of the spindle structure between two 
atoms but with negative eigenvalue region between them. 

We next examine the very weak Pd-H bond which is found in the four-fold bond in 
Pd9H22. There are two H atoms which have bonds between four Pd atoms and two of 
them are very weak (see Tables [l] and [2]) . Fig. [s] focuses on one of such H atom (H(17)) 
and shows one weaker Pd-H bond (Pd(3)-H(17)) and one stronger bond (Pd(4)-H(17)) 
(other two bonds, Pd(l)-H(17) and Pd(7)-H(17), are off this plane). The stronger one has 
the spindle structure and is similar to the one in the Pd-H-Pd bridging bond. Although 
the weaker one has a long bond distance (2.59 A) and negative eigenvalue region around 
the Lagrange point, there is a flow which connects Pd(3) and H(17) so we regard this as 
a bond. In other words, this bond is characterized by the pseudo-spindle structure. This 
classiflcation that shorter bonds are characterized by spindle structures and longer ones 
are by pseudo-spindle structures seems to hold rather in general for the hydrogenated Pd 
clusters we have investigated. In Sec. |3.H we have deflned two groups "A" and "B" for the 
Pd-H bonds where the former (the latter) has bond length shorter (longer) than about 
L9 A. In Sec. 3T, we have pointed out that they are on slopes with different inclination in 
the bond order v.s. bond length plot. From the viewpoint of the electronic stress tensor, 
group A is characterized by a spindle structure and the group B by a pseudo-spindle 
structure. 

Finally, we discuss the H-H bonds in PdsHio, Pd6Hi4 and PdyHie. Their bond lengths 
are respectively L56 A, L63 A and 2.03 A, which are longer than that of the H2 molecule. 
From the structural point of view, the H-H in PdsHio is located outside the Pd cage. The 
H-H in Pd6Hi4 is completely contained in the Pd cage and the one in PdyHie is marginally 
within the Pd cage. As before, we show the electronic stress tensor and tension for these 
bonds in Figs. [O} [TOl and [TTj In these flgures, there are more than three labelled atoms but 
they are all on the same plane (the distance between the plane and the labelled atoms are 



less than 0.01 A). In Fig. 11, left panel, there is a region with negative eigenvalue (shown 



by a blue circular region) at the lower-center part of the panel. This is caused by the 
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existence of Pd(6) but, since it locates slightly above this plane (about 0.4 A), it is not 
labelled. 

In these figures, we see that all the regions between Pd and H have a spindle structures, 
suggesting covalency of the Pd-H bonds. As for the H-H bonds, there is a spindle structure 



for the one in PdsHio (Fig. |9]). The H-H bonds in FdeRu and PdrHie (Figs. [TO] and [IT]) 

are characterized by the negative eigenvalue region and eigenvectors connecting H atoms, 
namely a pseudo-spindle structure. 

3.3 Improving bond order definition by surface integral 

In this section, we propose a way to improve our definition of bond order and apply it to 



the hydrogenated Pd clusters. As is described in Sec. |2.2[ our bond order is defined 
from the energy density evaluated at the Lagrange point. In the previous paper |12|, using 
many kinds of hydrocarbon molecules, it has been shown that this definition manifests a 
nice feature as a bond order. 

However, it is not difficult to imagine a type of chemical bond which cannot be well 
characterized by a single point between two atoms. This would be true for chemical bonds 
where spatially extended rf-orbitals are involved. Then, we consider it is worthwhile to 
investigate this issue by expanding our bond order definition by taking the surface integral 
of energy density instead of the energy density at the Lagrange point. 

Then, we need to determine the surface over which we integrate the energy density. 
The most natural choice would be a "Lagrange surface" |10| which is constructed from a 
family of lines which going out from a Lagrange point (if a Lagrange surface includes a 
Lagrange point). Namely, we define bond order of the bond between atoms A and B as 

(fcre^id) 

where Sxb denotes the Lagrange surface between atoms A and B. As is the cases of 
(Eq. ([5])) and 6^ (Eq. ([6])), we normalize by the value of hydrogen molecule. 

Unfortunately, however, this Lagrange surface is not so easy to define numerically. 
Hence we instead take the surface integral over the plane which includes a Lagrange point 
and is perpendicular to the axis connecting two atoms. Note that such a plane coincides 
with a Lagrange surface in the case of homonuclear diatomic molecules. 



The results of the h^^s) calculation are shown in Figs. 12 and 13, We see that h^^s) 
is calculated to be larger than for most of the bonds in the hydrogenated Pd clusters. 
This comes from the fact that the energy density distribution in those clusters is spatially 
extended relative to that of the hydrogen molecule. The ratio of h^(s) to is especially 



large for Pd^Pd bonds as shown in Fig. 13, This in turn is considered to be due to the 



d-orbitals of Pd atoms participating in the bonds. 
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4 Summary 



In this paper, we have apphed recently developed method to analyze electronic structure 
via electronic stress tensor to the small hydrogenated Pd clusters whose structures have 
been known (l][2]. From the results of ab initio electronic structure calculation of these 
clusters, we calculated quantities which are defined at each point in space: stress tensor, 
tension, energy density and regional chemical potential. The chemical bond is character- 
ized by the Lagrange point, where the tension vanishes and bond orders are defined by the 
energy density or chemical potential at that point. We have confirmed that thus defined 
bond orders are useful to classify the chemical bonds in Pd clusters as had been done for 
Pt clusters [l3j and AI4 clusters |17|. 

For some of the bonds, we have done more detailed analysis by drawing the eigenvalue 
of stress tensor, corresponding eigenvectors and tension vectors. In particular, the weak 
bonds with long bond distances suggested by the Lagrange point search in our bond order 
analysis are confirmed to have a fiow of eigenvectors connecting bonded atoms. They 
include four- fold Pd^H bonds and non-trivial H-H bonds. As for the H-H bond, we found 
that the one in PdsHio has a positive eigenvalue region, "spindle structure", showing 
covalency like that of the H2 molecule whereas those in Pd6Hi4 and PdyHie have negative 
eigenvalue regions, "pseudo-spindle structure". As for the Pd-H bond, it has been shown 
to be classified into two groups. One with shorter bond lengths which is characterized by a 
spindle structure and another with longer bond lengths characterized by a pseudo-spindle 
structure. 

We proposed a terminology "pseudo-spindle structure" in this paper but actually, there 
already has been such a structure found in our previous study for C2H2 |9|. The negative 
eigenvalue of C2H2 is caused by the compressive stress nearby the C nuclei. In general, 
the stress tensor has a large negative eigenvalue in radial direction in neighborhood of a 
nucleus due the dominance of the attractive Coulomb force. In the case of C2H2, the bond 
length is too short that the internuclear region is immersed under the atomic compressive 
stress [9j. This is a pseudo-spindle structure in a strong bond. We can say that we have 
found in this paper two more types of pseudo spindle structure. One is the pseudo-spindle 
structure associated with very long and weak H-H bond and Pd-H bond. Another is the 
pseudo-spindle structure associated with a bond between metallic atoms, Pd~Pd bond. 

We also have introduced an extension to our bond order. The modified definition uses 
the integration of the energy density over the "Lagrange surface" instead of the energy 
density at the Lagrange point. This modification makes the bond order greater, which 
reflects the contribution of spatially extended rf-orbitals to the bonds in the clusters. 

We believe that this study has provided another useful example of our stress tensor 
approach to chemical bonds. Further applications to other compounds, especially to those 
including transition metals, will solidify the basis of the stress tensor analysis and will 
deepen our understanding of chemical bonds. 
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Table 1: Number of bonding type in the hydrogenated Pd clusters. Pd-H (terminal), two- 
fold Pd H (bridging), three- fold Pd-H and four- fold Pd-H respectively count the number 
of H atoms who have one, two, three and four bonds between Pd atom. 



Cluster 


Pd-Pd 


Pd H 


Two-fold 


Three-fold 


Four-fold 


H-H 






(terminal) 


Pd-H 


Pd-H 


Pd-H 




Pd2H2 


1 





2 











PdsHs 


3 





1 


1 








Pd4H8 


3 


2 


6 











PdsHio 


4 





10 








1 


Pd6Hi4 


6 





12 


2 





1 


PdrHie 


12 


2 


13 


1 





1 


PdgHie 


8 





16 











Pd9H22 


13 


6 


10 


4 


2 






Table 2: The bond length of the Pd-H bonds participating in three-fold and four- fold 
Pd-H. See Fig. [T]for the label of H atom. 



Cluster 


Bond center 


Bond lengths [A] 


PdsHa 


H(4) 


1.69, 1.87, 1.87 


PdeHu 


H(10) 


1.69, 1.69, 1.69 




H(ll) 


1.69, 1.69, 1.69 


PdyHie 


H(19) 


1.75, 2.06, 2.06 


Pd9H22 


H(13) 


1.69, 1.76, 2.17 




H(14) 


1.72, 1.97, 2.03 




H(17) 


1.68, 1.71, 2.39, 2.59 




H(18) 


1.73, 1.97, 2.03 




H(20) 


1.68, 1.72, 2.37, 2.57 




H(21) 


1.70, 1.78, 2.11 
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Table 3: Linear and exponential fits to the bond order v.s. bond length data. We compute 
fits to three groups of data points separately: shorter Pd-H bonds (group A), longer Pd-H 
bonds (group B) and Pd-Pd bonds. In the fitting functions, x stands for the bond length. 
Xred stands for the reduced of the fit. 



Bond Types 


be — ax + b 


be — cexp{—dx) 


bfj, — ax -\- b 


bfj, — cexp{—dx) 


Pd-H (A) 


a = -1.43 
6 = 2.81 
xL = 5.9 X 10-4 


c = 2.42 X 10^ 
d = 3.82 
xL = 1-3 X 10-4 


a = -1.40 
b = 3.33 
X?ed = 7.3 X 10-4 


c= 11.5 
d= 1.47 
xL = 8.3 X 10-4 


Pd-H (B) 


a = -0.173 
b = 0.487 

X?ed = 1-3 X 10-4 


c = 8.53 
d = 2.04 

X?ed = 3.6 X 10-5 


a = -0.270 
6= 1.22 

Xla = 2.4 X 10-4 


c= 1.66 
d = 0.448 

Xla = 2.3 X 10-4 


Pd-Pd 


a = -0.346 
b = 1.08 
X?ed = 3.7 X 10-5 


c = 1.40 X 10^ 
d = 2.54 
xL = 2.6 X 10-5 


a = -0.931 
b = 3.36 
X?ed = 4.8 X 10-4 


c= 17.9 
d^ 1.13 
= 6.0 X 10-4 
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Figure 1: Optimized structures and bond order for hydrogenated Pd clusters: (a) Pd2H2, 
(b) Pd3H2, (c) Pd4H8, (d) PdsHio, (e) Pd.Ru. (f) PdrHig, (g) PdgHie and (h) Pd9H22. 
The bonds are drawn at which Lagrange points are found and our energy density based 
bond order (eq. ^) is shown by color. 
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Figure 3: Comparison of the relation between bond length and various bond orders: the 
Wiberg bond index (blue square) , atom- atom overlap- weighted N AO bond order (magenta 
circle), the Mayer's bond order (light-blue triangle) and our bond orders (red down- 
triangle) and 6^ (green diamond). 



16 



Pd-H(terminal) 
Pd-H (two-fold) 
Pd-H(three-fold) 
Pd-H(four-fold) 
Pd-Pd(PdJH 
Pd-PdM) 
H-H(PdnH^) 



2 2.2 2.4 2.6 2.8 

Bond length[A] 



Figure 4: The relation between bond length and energy density bond order in the bare 
Pd clusters and hydrogenated Pd clusters. We classify bonding types as terminal Pd-H 
(light-blue cross), two- fold Pd-H (green asterisk), three- fold Pd-H (blue square), four- fold 
Pd-H (black triangle), Pd-Pd in hydrogenated Pd clusters (magenta circle), Pd-Pd in 
bare Pd clusters (red diamond), and H-H (red fiUed-triangle) . Note that the H2 molecule 
would locate at (0.748, 1.000). 



□ □ 



Pd-H(terminal) 
Pd-H (two-fold) 
Pd-H(three-fold) 
Pd-H(four-fold) 
Pd-Pd(PdJH J 
Pd-Pd(Pd„) 
H-H(PdnHJ 




2 2.2 2.4 

Bond length[A] 



Figure 5: Same as Fig. lilbut with chemical potential bond order b 
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Figure 6: The H-H bond in H2. The largest eigenvalue of the stress tensor and correspond- 
ing eigenvector (left panel) and tension (right panel) are plotted on the plane including 
the bond axis. On the left panel, contours of 0.1 and —0.1 are shown by white dotted 
lines. On the right panel, the normalized tension vectors whose norm is expressed by the 
color of the arrows are shown and the Lagrange point is marked by a black diamond. 
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Figure 7: The Pd-H-Pd bridging bond in Pd6Hi4 is shown as the largest eigenvalue of the 
stress tensor and corresponding eigenvector (left panel) and tension (right panel) on the 
plane including the labelled atoms. See Figs, [l] or [2] for the number in the label. On the 
right panel, the Lagrange point is marked by a black diamond. 
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Figure 8: The four-fold Pd-H bond in Pd9H22 is shown in the similar manner as Fig. [7| 
Two of four Pd-H bonds from H(17) are shown on this plane. 
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Figure 9: The H-H bond in PdsHio is shown in the similar manner as Fig. [7} 
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Figure 10: The H H bond in Pd6Hi4 is shown in the similar manner as Fig. [7| 
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Figure 11: The H-H bond in PdrHig is shown in the similar manner as Fig. [7| 
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Figure 12: Comparison of the relation between bond length and bond orders. Bond order 
defined from the energy density at the Lagrange points (6^) and one defined from surface 
integral of the energy density (b^^s))- 



4.5 



3.5 - 



g 2.5 

CO 

2 



Pd-H(terminal) 
Pd-H (two-fold) 
Pd-H(three-fold) 
Pd-H(four-fold) 
Pd-Pd(Pd H J 
Pd-Pd(Pclj 
H-H(Pd,HJ 



□ * 



□ 

5g 




1.5 



0.5 



1.4 



1.6 



1.8 



2.2 2.4 

Bond length[A] 



2.6 



2.8 



Figure 13: The ratio of b^^s) to 6^. 
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